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Freezing of charge-stabilized colloidal suspensions and rel- 
ative stabilities of crystals and quasicrystals are studied us- 
ing thermodynamic perturbation theory. Macroion interac- 
tions are modelled by effective pair potentials combining elec- 
trostatic repulsion with polymer-depletion or van der Waals 
attraction. Comparing free energies - counterion terms in- 
cluded - for elementary crystals and rational approximants to 
icosahedral quasicrystals, parameters are identified for which 
one- component quasicrystals are stabilized by a compromise 
between packing entropy and cohesive energy. 

PACS numbers: 82.70.Dd, 83.70.Hq, 61.44.Br, 05.70.Fh 



Suspensions of mesoscopic-sized (1 — 1000 nm) colloidal 
particles dispersed in a fluid medium are of both prac- 
tical and fundamental interest Q. Beyond traditional 
relevance to the chemical, food, and pharmaceutical in- 
dustries, the remarkable mechanical, thermal, and optical 
properties of these typically soft materials raise intriguing 
prospects for novel applications, such as optical switch- 
ing devices Recent scientific interest has been driven 
largely by advances in sample preparation, scattering, 
and imaging techniques As well-characterized mod- 
els of atomic systems, colloidal suspensions offer valuable 
insight into the basic link between microscopic interpar- 
ticle interactions and macroscopic phase behavior in con- 
densed matter. More so than in atomic systems, colloidal 
interactions are eminently tunable through experimental 
control of particle size, charge, and surface chemistry, 
as well as properties of the suspending medium, such as 
polarizability and salt concentration. A correspondingly 
rich variety of thermodynamic phases has been observed, 
including the vapor, liquid, crystal, and glass phases com- 
mon to atomic systems. A notable exception is the qua- 
sicrystal phase, which to date is known to exist only in 
certain binary and ternary metallic alloys [Q. The pur- 
pose of this Letter is to demonstrate, using general sta- 
tistical mechanical methods, the possibility of thermody- 
namically stable one- component colloidal quasicrystals. 
If experimentally observed, such systems would represent 
a fundamentally new type of colloidal order and permit 
the first real-space imaging of quasicrystals. 

We focus specifically on suspensions of charged spher- 
ical colloidal macroions and oppositely charged counte- 
rions dispersed in a deionized (salt-free) solvent. The 
macroions are modelled as hard spheres of diameter a 
and uniform surface charge Ze, the counterions as point 



particles of elementary charge e, and the solvent as a 
continuous medium of dielectric constant e. Charge neu- 
trality relates the average macroion and counterion num- 
ber densities via Zpm — Pc- Electrostatic interactions 
between charged macroions in suspension are commonly 
modelled by an effective hard-sphere/screened-Coulomb 
pair potential: 
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where r is the separation between macroion centers and 
K = -y/ At: Z pra I eksT is the Debye screening constant. 
Equation (^ was first derived in the dilute regime by Der- 
jaguin. Landau, Verwey, and Overbeek The screened- 
Coulomb form has since been shown to be accurate also 
for concentrated suspensions, with Z replaced by an ef- 
fective charge [^. At finite macroion volume fractions 
rj — [t: / Q) pmO^ , the volume available to the counterions 
is smaller than the total volume F by a factor (1 — rj). 
As a result, the effective counterion density is increased 
by a factor 1/(1 — 77), enhancing screening Q according 
io k — -y/ Att Z pm / eksT {1 — rj). Henceforth, we assume 
an electrostatic potential of the form of Eq. (|^), but with 
this modified screening constant. 

Aside from electrostatic repulsions, interactions be- 
tween macroions also can include an attractive compo- 
nent. Attractions may arise, for example, in the presence 
of free (nonadsorbing) polymer via a polymer-depletion 
mechanism . For macroion surface-surface separations 
smaller than the characteristic polymer coil diameter, de- 
pletion of polymer from the intervening space creates an 
osmotic pressure imbalance drawing the macroions to- 
gether. The attractive energy is directly proportional to 
the polymer osmotic pressure lip and the mutual over- 
lap of excluded volume shells. For coils of sufficiently 
small radius of gyration Rg relative to cr, the effect is 
described by a pair potential that is nonzero only in the 
range a < r < a 



2Rg , where it has the form 
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with X = r/a, ^ = 2Rg/a, and c = {T:/6)Ilpa\l + • 
For size ratios ^ > 0.16, triplet overlaps entail three-body 
interactions. However, for ^ ~ 0.25 (our case below) 
these are small enough to be ignored M. Attractions 
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also may arise from fluctuating dipole-dipole forces, as 
described by the London- van der Waals pair potential H] 
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where the Hamaker constant A is related to the macroion 
and solvent polarizabilities (refractive indices). The neg- 
ative divergence at contact (r = a) can cause irreversible 
coagulation unless the macroions are stabilized, e.g., ster- 
ically by adsorbing or grafting a layer of polymer to their 
surfaces. For simplicity, we model this by a cut-off at a 
distance r = Rc, roughly 5 — 10 % larger than a. 

Figure 1 depicts effective macroion-macroion pair po- 
tentials that combine electrostatic repulsion with each of 
the two types of attraction discussed above. The length 
scales of the interactions are set by the macroion core 
diameter, the polymer coil radius (or van der Waals cut- 
off distance), and the Bjerrum length Xb = e'^/eksT (at 
temperature T). These pair potentials are now taken as 
input to a statistical mechanical theory for the Helmholtz 
free energy. The macroion contribution to the free energy 
is determined by means of thermodynamic perturbation 
theory |ll[] . (An earlier such study by Gast et al. [|l2| con- 
sidered freezing only into an FCC crystal.) The theory 
proceeds by splitting the pair potential into short-range 
reference and longer-range perturbation potentials. The 
natural separation here is into a hard-sphere (HS) refer- 
ence potential (j)HS (r) and a perturbation potential (f>p (r) 
that combines an attractive well with a repulsive barrier 
and long-range tail. To first order in 0p(r), the macroion 
free energy separates correspondingly into reference and 
perturbation terms: 

Fm[Pmir)] = FHs[Prn{r)] 

poo 

+ 2np,nNm / dr' r'^ gHs{r'\ [Prn{v)]) 0p(r'), (4) 
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Nm = PmV being the macroion number, and -Ff/s[Pm(r)] 
and gnsif", [Pm(r)]) the free energy and radial distribu- 
tion function (RDF), respectively, of the HS reference 
system, both dependent on the equilibrium macroion den- 
sity distribution p„i{v). Accuracy of first-order perturba- 
tion theory is assured if fluctuations in the total pertur- 
bation energy remain small relative to ksT [ pT| , a con- 
dition generally well satisfied here, where variations in 
(/)(r) are small for distances over which the RDF varies 
appreciably. 

For the fluid phase, with spatially constant macroion 
density, the macroion free energy is calculated via the 
uniform limit [jOm(r) pm] of Eq. (^), using the essen- 
tially exact Carnahan-Starling and Verlet-Weis forms |l^] 
for FHsipm) and gHsir; Pm), respectively. For the solid 
phase, classical density-functional methods are ap- 
plied. The reference free energy is obtained by minimiz- 
ing a, functional J- H s[Pm{i^)] with respect to the macroion 



density distribution pm(r), or in practice the width pa- 
rameter a in the Gaussian parametrization 
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the sum running over the lattice sites R of a spec- 
ified solid structure. For non-overlapping Gaus- 
sians, the configurational entropy (ideal-gas) part of 
the functional is given exactly by J^id/NmkBT — 
{i/2)\n{aa'^ /it) — 5/2. The excess part is determined 
by means of the modified weighted-density approxima- 
tion (MWDA) |]l^, which reasonably describes HS solids. 
This maps the excess part of J-Hs[Pm{T^)\l onto its 
fluid counterpart, evaluated at an effective density, via 
•^il'^^'^KWl/^m = /e,(/5), where/,, is the fluid ex- 
cess free energy per particle and p is a weighted average 
of Pm (r) that incorporates exact two-particle correlations 
and a subset of higher-order correlations. The reference 
RDF - an angular average of the two-particle density 
- is calculated using an approach recently proposed by 
Rascon et al. . This corrects the first peak for nearest- 
neighbor correlations - by fixing the contact value (via 
the virial theorem), coordination number, and first mo- 
ment - and treats higher-order peaks in a mean-field fash- 
ion. The approximation compares excellently with sim- 
ulation data for HS crystals and has been successfully 
applied to Lennard- Jones and square- well solids . 
The counterion contribution to the free energy is [lO| ] 
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comprising the entropy of the counterions and the in- 
teraction of a macroion with its own counterion screen- 
ing cloud. Although independent of macroion structure, 
the counterion free energy is manifestly state-dependent 
and, under conditions of high charge and low salt, can 
profoundly influence phase behavior. 

Applying the above procedure, we have calculated to- 
tal free energies for a variety of elementary crystals and 
model quasicrystals. The crystalline structures exam- 
ined include face-centered and body-centered cubic (FCC 
and BCC), tetragonal, orthorhombic, and rhombohedral 
Bravais lattices, and the hexagonal close-packed (HCP) 
structure. The model quasicrystalline structures are ra- 
tional approximants (RA) to icosahedral quasicrystals 
generated by projecting a six-dimensional hypercubic lat- 
tice onto the three-dimensional physical space and ap- 
proximating, in the perpendicular space, the golden mean 
T by a rational number r„ = Fn+i/ Fn, where Fn is a term 
in the Fibonacci sequence. This procedure, together with 
a one-component occupation of the lattice sites that opti- 
mizes packing efficiency , yields periodic lattices with 
large unit cells and local order closely approximating that 
of aperiodic quasicrystals. The first four RAs, denoted 
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1/1, 2/1, 3/2, and 5/3, have unit cells containing, respec- 
tively, 20, 108, 452, and 1904 atoms and maximum HS 
volume fractions of 0.5020, 0.6400, 0.6323, and 0.6287 
(c/. 0.6288 in the quasiperiodic limit n ^ oo and 0.7405 
for FCC and HCP crystals). The RDFs illustrated in 
Fig. 2 show that nearest-neighbor distances for the RAs 
are about 5 % shorter than for the crystals. For compu- 
tational ease we focus on the 2/1 structure, but check to 
ensure qualitatively similar results for the higher-order 
approximants. 

From previous work p^ , packing efficiency is known 
to dominate the stability of HS solids, the RAs being al- 
ways metastable (higher in free energy) relative to the 
crystals. On the other hand, because of their shorter 
nearest-neighbor distances, the RAs can gain energeti- 
cally from short-range attractive interactions. In fact, 
for potential wells of appropriate depth (a few fc^T) and 
width ( 10 — 30 % of (t), their free energies can be lower 
than those of simple crystals. As shown above, such at- 
tractions can be produced in colloidal systems by adding 
polymer of appropriate size and concentration [Eq. (H)] 
or by tuning refractive indices [Eq. (^)]. 

Assuming a room-temperature, salt-free, aqueous sol- 
vent (Ab — 0.72 nm), we have surveyed macroion param- 
eters, a and Z, and parameters defining attractive in- 
teractions - ^ and Hp for index-matched colloid-polymer 
mixtures, Rc and A for non-index-matched suspensions - 
computed free energies for each solid structure, and iden- 
tified parameter combinations tending to favor quasicrys- 
tal stability. For illustration, Fig. 3 displays the free ener- 
gies of a non-index-matched suspension characterized by 
cr = 50 nm, Z = 150e, Rc/a = 1.058, and A/keT = 25. 
Hamaker constants of this magnitude, though high for 
polymeric colloids, are typical for metallic colloids, such 
as Ag, Au, and Cu. Furthermore, charges of this mag- 
nitude correspond to counterion concentrations of order 
10"'^ mol/1, far exceeding the background concentration 
of 10~^ mol/1 for H+ and OH^ ions from the dissociation 
of water. Relative stabilities of the fluid and compet- 
ing solid structures are assessed by constructing Maxwell 
common tangents, ensuring equality of chemical poten- 
tials and pressures in coexisting phases. Remarkably, 
over a range of densities intermediate between dense fluid 
and close-packed crystals, the quasicrystalline structures 
are predicted to be thermodynamically stable. The sta- 
bility can be traced to a competition between two fac- 
tors, namely packing efficiency, favoring the FCC/HCP 
crystal, and nearest-neighbor coordination, favoring qua- 
sicrystalline order. This is best seen by comparing Figs. 1 
and 2 and observing that the attractive well of (p{r) is 
more commensurate with the first peak of gusi''^) for the 
2/1 RA than for the FCC crystal. Stable quasicrystals 
thus emerge as a compromise between packing entropy 
and cohesive energy in salt-free colloidal suspensions with 
a concentrated counterion background. 

Applying the same approach to colloid-polymer mix- 



tures, quasicrystalline stability is again predicted for suf- 
ficiently high polymer osmotic pressure. Repeating the 
coexistence analysis over a range of osmotic pressures, we 
have mapped out phase diagrams in the Ilp-pm plane. A 
typical example is presented in Fig. 4. At lower osmotic 
pressures, where the interactions are purely repulsive, the 
fluid freezes into a close-packed crystal, FCC and HCP 
being essentially degenerate. Beyond a threshold osmotic 
pressure, however, freezing occurs into the quasicrystal, 
which remains stable over a significant range of densi- 
ties. At the threshold pressure (here Hpa^ /ksT ~ 220), 
which we estimate to be near cross-over between dilute 
and semi-dilute polymer regimes, the fluid, crystal, and 
quasicrystal coexist at a triple point. At higher densities, 
packing efficiency prevails and the quasicrystal makes 
a structural transition to the FCC/HCP crystal. The 
small relative density differences justify our assumption 
of equal polymer concentrations in coexisting phases. Fi- 
nally, we emphasize that the predicted region of qua- 
sicrystal stability is robust with respect to variation of 
parameters, widening and shifting to higher Hp with in- 
creasing a and Z, and that its existence depends not on 
specific details, but only the qualitative form, of the in- 
teractions. 

Summarizing, using thermodynamic perturbation the- 
ory, we have investigated freezing transitions and as- 
sessed the relative stabilities of crystalline and qua- 
sicrystalline solids in charge-stabilized colloidal suspen- 
sions. For effective macroion interactions combining 
electrostatic repulsion with either polymer-depletion or 
van der Waals attraction, system parameters and ther- 
modynamic states have been identified for which one- 
component icosahedral quasicrystals are predicted to be 
thermodynamically stable over a significant range of den- 
sities. Such systems and interactions should be experi- 
mentally achievable, raising prospects for observation of 
stable colloidal quasicrystals. 

We thank Prof. J. Hafner and Dr. A.M. Denton for 
helpful discussions, and Drs. M. Krajci and M. Windisch 
for supplying the rational approximant structures. 
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FIG. 1. Effective pair potentials for charge-stabilized colloidal suspensions with a — 50 run, Z = 150e. Solid curve: 
index-matched suspension mixed with polymer of Rg/cr = 0.125, Upa^ /ksT = 250. Dashed curve: non-index-matched with 
Rc/a = 1.058, A/UbT = 25. 
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FIG. 3. Helmholtz free energy per unit volume vs. macroion density (in reduced units) for a colloidal suspension interacting 
via pair potential of Fig. 1 (dashed curve) . 




FIG. 4. Phase diagram of polymer osmotic pressure vs. macroion density for a colloid-polymer mixture interacting via pair 
potential of Fig. 1 (solid curve). Horizontal tie lines connect corresponding points on coexistence curves. 
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